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Abstract 

Over the last decades strong efforts have been made 
to apply new spaceborne technologies to the study 
and possible forecast of strong earthquakes. In this 
study we use ASTER/TERRA multispectral satellite 
images for detection and analysis of changes in the 
system of lineaments previous to a strong earthquake. 
A lineament is a straight or a somewhat curved fea- 
ture in an image, which it is possible to detect by 
a special processing of images based on directional 
filtering and or Hough transform. "The Lineament 
Extraction and Stripes Statistic Analysis" (LESSA) 
software package, developed by Zlatopolsky (1992, 
1997). We assume that the lineaments allow to de- 
tect, at least partially, the presence ruptures in the 
Earths crust, and therefore enable one to follow the 
changes in the system of faults and fractures associ- 
ated with strong earthquakes. We analysed 6 earth- 
quakes occurred in the Pacific coast of the South 
America and one earthquake in Tibet, Xizang, China 
with the Richter scale magnitude > 5.2 Mw. They 
were located in the regions with small seasonal vari- 
ations and limited vegetation to facilitate the track- 
ing of features associated with the seismic activity 
only. It was found that the number and orientation 
of lineaments changed significantly about one month 
before an earthquake approximately, after that the 
system gradually returns its initial state. This ef- 
fect increases with the earthquake magnitude, and 



it is much more easily detectable in case of conver- 
gent plate boundaries (for example, Nazca and South 
American plates) . The results obtained open a possi- 
bility to develop a methodology able to evaluate the 
seismic risk in the regions with similar geological con- 
ditions. 

1 Introduction 

Throughout the world, devastating earthquakes con- 
stantly occur with little or no advance warning. 
Brune (1979) proposed that earthquakes may be in- 
herently unpredictable since large earthquakes start 
as smaller earthquakes, which in turn start as smaller 
earthquakes, and so on. In his model, most of the 
fault is in a state of stress below that required to ini- 
tiate slip, but it can be triggered and caused to slip 
by nearby earthquakes or propagating ruptures. Any 
precursory phenomena will only occur when stresses 
are close to the yield stress. However, since even 
small earthquakes are initiated by still smaller earth- 
quakes, in the limit, the region of rupture initiation 
where precursory phenomena might be expected is 
vanishingly small. Even if every small earthquake 
could be predicted, one still faces the impossible task 
of deciding which of the thousands of small events 
will lead to a runaway cascade of rupture composing 
a large event. 

Nevertheless, the discussion about the possibility 
of earthquake forecast continues to be open, and a 



wide spectrum of new spaceborne technologies for 
the earthquake study and forecast appeared during 
the last decades. The main advantage of spaceborne 
technologies is the ability to cover big territories and 
areas with difficult access. The list of these tech- 
nologies is very large. As an example it is possi- 
ble to mention the measurements of different iono- 
spheric precursors of earthquakes including changes 
in electromagnetic ELF radiation (Serebryakova et 
al., 1992, Gokhberg et al., 1995), and ionospheric 
electron temperature (Sharma et al., 2006), and den- 
sity (Trigunait et al., 2004) (see Pulinets at al., (2003) 
for a review). Many efforts have been concentrated in 
the study of the ground deformation using the satel- 
lite radar interferometry, that makes it possible to 
determine the location and amount of coseismic sur- 
face displacements (see for example Satybala, 2006; 
Schmidt and Bergmann, 2006, Lasserre et al., 2005, 
Funning et al. , 2005) . The IR satellite thermal imag- 
ing data were used to study pre-earthquake ther- 
mal anomalies (Ouzounov and Freund, 2004). The 
anomalies in the surface latent heat flux data were 
also detected a few days prior to coastal earthquakes 
(Cervone et al., 2005, Singh and Ouzounov, 2003; 
Dey at al., 2004). 

During last years, significant progress has been 
reached in the understanding how the complex set of 
phenomena, related to the earthquake gestation is re- 
flected, at least partially, in the geological lineaments. 
In particular, Cotilla-Rodriguez and Cordoba-Barba 
(2004) studied the morphotcctonic structure of the 
Iberian Peninsula and showed that the main seismic 
activity is concentrated on the first- and second rank 
lineaments, and some of important epicenters are lo- 
cated near the lineament intersections. Stich et al., 
(2001) found from the analysis of 721 earthquakes 
with magnitude between 1.5 and 5.0, that the epi- 
centers draw well-defined lineaments and show two 
dominant strike directions N120-130E and N60-70E, 
which are coincident with the known fault system of 
the area. Distances within multiplcts (typically sev- 
eral tens of meters) are smaller than the fracture radii 
of these events. Carver et al. (2003) have used the 
SRTM and Landsat-7 digital data and paleoseismic 
techniques to identify active faults and evaluate seis- 
mic hazards on the northeast coast of Kodiak Island, 



Alaska. 

Arellano et al. (2004ab, 2005) studied the changes 
in the lineament structure caused by a 5.2 Richter 
scale magnitude earthquake occurred January 27, 
2004 in southern Peru. During last years this re- 
gion is studied intensively using the ground based 
seismic network (Comte et al., 2003; David et al., 
2004; Legrand, 2005) as well as GPS and SAR 
interferometry data (Campos et al., 2005). The 
ASTER/TERRA high resolution multispectral im- 
ages 128 and 48 days before and 73 days after the 
earthquake were used. It was shown that the linea- 
ment system is very dynamical, and significant num- 
bers of lineaments appeared between four and one 
month before the earthquake. They also studied the 
changes in stripe density fields. These fields represent 
the density of stripes, calculated for each direction 
as a convolution between the corresponding circular 
masks and the image. The stripe density field resid- 
uals showed the reorientation of stripes, which agrees 
with the dilatancy models of earthquakes. These fea- 
tures disappear in the image obtained two months af- 
ter the earthquake. Analysis of the similar reference 
area, situated at 200 km from the epicenter, showed 
that in the absence of earthquakes both lineaments 
and stripe density fields remain unchanged. Simi- 
lar results were obtained later by Bondur and Zverev 
(2005) due to analysis of MODIS (TERRA) images 
of earthquake in California. 

Singh V.P and R.P. Singh (2005) used the lin- 
eament analysis to study changes in stress pattern 
around the epicenter of Mw=7.6 Bhuj earthquake. 
This earthquake occurred 26 January 2001 in In- 
dia. Indian Remote Sensing (IRS-1D) LISS data 
were used. The lineaments were extracted using high 
pass filter (Sobel filter in all directions). The results 
obtained also confirm that the lineaments retrieved 
from the images 22 days before the earthquake dif- 
fer from the lineaments obtained 3 days after the 
earthquake. It was assumed that they are related 
to fractures and faults and their orientation and den- 
sity give an idea about the fracture pattern of rocks. 
The results also show the high level of correlation 
between the continued horizontal maximum compres- 
sive stress deduced from the lineament and the earth- 
quake focal mechanism. 
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Studies of lineament dynamics can also contribute 
to better understanding of the nature of earthquakes. 
To date significant number of theories has been de- 
veloped to explain how an earthquake occur. One 
of the oldest is the elastic rebound theory, proposed 
by Harry Reid after the California 1906 earthquake 
(Reid, 1910). It is based on the assumption that 
the rocks under stress deform elastically, analogous 
to a rubber band. Strain builds up until either the 
rock break creating a new fault or movement oc- 
curs on an existing fault. As stored strain is re- 
leased during an earthquake, the deformed rocks "re- 
bound" to their undeformed shapes. The magnitude 
of the earthquake reflects how much strain was re- 
leased. The seismic gap hypothesis states that strong 
earthquakes are unlikely in regions where weak earth- 
quakes are common and the longer the quiescent pe- 
riod between earthquakes, the stronger the earth- 
quake will be when it finally occur (see Kagan and 
Jackson, 1995, and references therein). The compli- 
cation is that the boundaries between crustal plates 
are often fractured into a vast network of minor faults 
that intersect the major fault lines. When an earth- 
quake relieves the stress in any of these faults, it may 
pile additional stress on another fault in the network. 
This contradicts the seismic gap theory because a se- 
ries of small earthquakes in an area can then increase 
the probability that a large earthquake will follow. 

The theory of dilatancy states that an earthquake 
develops similarly to the rupture of a solid body 
(Whitcomb ct al., 1973; Scholz et al., 1973; Griggs et 
al., 1975). This approach has a physical basis in lab- 
oratory studies of rock samples, which showed that 
when rocks are compressed until they fracture, a di- 
latancy often occurs for a short time interval imme- 
diately before failure (Scholz, 1968). Mjachkin et al. 
(1975ab) modified the dilatancy approach and formu- 
lated the theory of unstable avalanche crack forma- 
tion. The model is based on the two phenomena: in- 
teraction between the stress fields of the cracks, and 
the localization of the process of the crack forma- 
tion. The number and size of cracks increases grad- 
ually under the action of tensions below a critical 
value. When the density of cracks reaches some crit- 
ical value, the rock breaks very quickly. This pro- 
cess develops due to merging of cracks as a result 



of interaction between their stress fields. However, 
the larger cracks have more probability to interact, 
and it supposes that a small number of large cracks 
is gradually formed, and their merging leads to the 
macro-destruction. During the earthquake gestation, 
a gradual increase in number and size of cracks oc- 
cur in the whole volume of rock under compression. 
When the crack density reaches a critical value, the 
barriers between cracks are destroyed, and the ve- 
locity of deformation increases. Finally, an unstable 
deformation develops and localizes in a narrow zone 
of future macro-rupture, the cracks orient along the 
future macro rupture, and a macro-crack is formed, 
producing an earthquake. 

However, this model was modified recently by in- 
troducing a concept of self-organized criticality, pro- 
posed by Bak et al. (1988) for description of the be- 
havior of complex systems. Applied to earthquakes, 
this approach describes an interaction between the 
ruptures of different rank and the collective effects 
of rupture formation before a strong earthquake (for 
example Varnes, 1989; Keilis-Borok, 1990; Sammis 
and Sornette, 2002). A wide area around the future 
epicenter reaches a metastable state, and the system 
turns to be very sensitive to small external actions. 
The concept of SOC does not contradict the concept 
of dilatancy. However, it assumes that a significantly 
greater region is involved during the last stages of 
earthquake preparation than the dilatancy theories 
imply. 

Unfortunately, the main processes leading to an 
earthquake develop deep inside the crust, and there 
is no way to realize direct measurements of any quan- 
tity. The unique possibility we have is to search 
for traces of these processes disseminated over the 
Earth's surface. In this context, the lineament anal- 
ysis could convert in the future in one of power tools 
for earthquake study, complementing other ground- 
based and satellite studies. Nevertheless, despite 
promising results obtained, many important ques- 
tions continue to be present. It is necessary to un- 
derstand, whether the lineament system is always af- 
fected by earthquake? How early before an earth- 
quake is this alteration manifested? How is it related 
to the earthquake magnitude and depth? How dif- 
ferent is it in case of different kinds of plate borders? 
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This study represent a first step in the search of some 
answers. 

2 Instrumentation and Data 
Analysis 

For this study we used the the images from the Ad- 
vanced spaceborne Thermal Emission and Reflection 
Radiometer (ASTER) onboard the TERRA satel- 
lite. The satellite was launched to a circular solar- 
synchronous orbit with altitude of 705 km. The ra- 
diometer is composed by three instruments: Visible 
and Near Infrared Radiometer (VNIR) with 15 m res- 
olution (bands 1-3), Short Wave Infrared Radiometer 
(SWIR) with 30 m resolution (bands 4-9) and Ther- 
mal Infrared Radiometer TIR with 90 m resolution 
(bands 11-14) which measure the reflected and emit- 
ted radiation of the Earths surface covering the range 
0.56 to 11.3 fim (Abrams, 2000). 

The images were processed using the Lineament 
Extraction and Stripes Statistic Analysis (LESSA) 
software package (Zlatopolsky, 1992, 1997), which 
provides a statistical description of the position and 
orientation of short linear structures through detec- 
tion of small linear features (stripes) and calculation 
of descriptors that characterize the spatial distribu- 
tion of stripes. The program also makes it possible 
to extract the lineaments - straight lines crossing a 
significant part of the image. To make this extrac- 
tion, a set of very long and very narrow (a few pixels) 
windows (bands), crossing the entire image in differ- 
ent directions, was used. In each band the density 
of stripes, the direction of which is coincident with 
the direction of the band, is calculated. When the 
density of stripes overcomes a pre-established thresh- 
old, the chain of stripes along the band is considered 
as a lineament. The value of threshold depends on 
the brightness of the image, relief, etc. and is es- 
tablished empirically. Previous studies showed that 
lineaments, extracted from the image by applying the 
LESSA program, are strongly related to the main lin- 
eaments, obtained from the geomorphological studies 
(Zlatopolsky, 1992, 1997). The details about the ap- 
plication of LESSA package for earthquake studies is 



given in (Arellano et al., 2005). 

During this study we analysed 5 earthquakes, oc- 
curred in the in the Pacific coast of the South Amer- 
ica and one earthquake occured in Himalaya, China. 
Table 1 resumes main characteristics of these earth- 
quakes, indicating the date, country, geographic co- 
ordinates, magnitude, and depth of the earthquake. 
Also the the ASTER images available for each earth- 
quake are indicated, for example -126 means that the 
image 126 days before earthquake was used. The last 
column indicates that in all South American earth- 
quakes number and orientation of lineaments suffered 
changes before the earthquake. In case of China 
earthquake, we can not give a defnite answer, be- 
cause unfortunately the key images tens day before 
the earthquake were covered by clouds in appoxi- 
mately 50%, that made the lineament analysis diffi- 
cult (last two lines, two areas cvering the hipocenter 
and close to hipocenter). Neverthless, more sophys- 
ticated technique based on analysis of stipe density 
fields was able to detect the alterations in these fields 
related to the earthquake. The methodology of this 
analysis is given in (Arellano et al., 2006). Currently 
we are preparring a manuscript dedicated especially 
to the analysis of this event. 

To illustrat the results obtained we give as an 
example a detailed analysis of 7.8 Mw earthquake, 
which took place June 13, 2005 in northern Chile 
close to Arica (see Figure 1). The hipocenter was sit- 
uated at 115 km deep in the crust. The coordinates 
were -19.99o LAT, -69.197o LONG. In the top, a se- 
ries of four band 3 ASTER (VNIR) images around the 
hipocenter area are shown. It is possible to see, that 
the presence of clouds was low. The second line con- 
tains the images showing the systems of lineaments, 
obtainend from the images above using the LESSA 
programm with a threhhold 120 (Zlatopolsky, 1992, 
1997). It is posible to see clear time evolution of lin- 
eaments, experimenting strong increase in the num- 
ber of lineaments 5 days before the earthquake. The 
third and fourth lines quanify this effect by calculaing 
the rose-diagramms and the Radon transforms. Re- 
orientation of lineaments can be taken as an indirect 
evidence in favour to the theory of dilatney. Never- 
theless, it is necessary to make more detailed studies 
to make definitive conclusions. 
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Place 


Date 


Magn.Mw. 


Depth 


Lat. 


Long. 


Images available 


Changes in lineaments 


Chile 


6/13/05 


7.8 


115.6 


-19.99 


-69.197 


-126, -69, -5, +139 


yes 


Chile 


9/17/03 


5.8 


127.1 


-21.467 


-68.325 


-138, -12, +86 


yes 


Peru 


10/17/05 


5.8 


123 


-17.775 


-69.486 


-51, -3, +4, +132 


yes 


Chile 


6/17/04 


5.7 


115.4 


-21.246 


-68.372 


-188, -35, -28, -19, +100 


yes 


Peru 


01/27/04 


5.2 


120.1 


-17.69 


-70.65 


-128,-48,+73 


yes 


China 


4/19/06 


5.7 


33.1 


31.6 


90.4 


-132, -116, -52, -20*, -4* 


? 



Table 1 : Main characteristics of earthquakes analyzed 
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February 7, 2005 April 5, 2005 June 8, 2005 October JO, 2(105 

(-126 days) (-69 days) (-5 day*) (+B9 days) 




Figure 1: From top to bottom: ASTER band 3 images around the June 13, 2005 earthquake (Chile, Arica). 
Systems of lineaments extracted from these images, using LESSA sofrware. 
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3 Discussion and conclusions 

In this study we used the multispectral satellite im- 
ages from ASTER/TERRA satellite for detection 
and analysis of lineaments in the areas around strong 
earthquakes with magnitude more than 5 Mw. A lin- 
eament is a straight or a somewhat curved feature in 
an image, which can be detected by a special process- 
ing of images, based on directional filtering and/or 
Hough transform. It was established that the sys- 
tems of lineaments are very dynamical. By analyzing 
5 events of strong earthquakes, it was found that a 
significant number of lineaments appeared approx- 
imately one month before an earthquake, and one 
month after the earthquake the lineament configura- 
tion returned to its initial state. These features were 
not observed in the test areas, situated hundreds kilo- 
meters away from the earthquake epicenters. 

The main question is how the lineaments extracted 
from images of 15-30 m (ASTER) in resolution are 
able to reflect the accumulation of stress deep in the 
crust given that the ground deformations associated 
with these phenomena are about a few centimeters? 
The nature of lineaments is related to the presence of 
faults and dislocations in the crust, situated at dif- 
ferent depth. If a dislocation is situated close to the 
surface, the fault appears as a clear singular linea- 
ment. In the case of a deep located fault, we ob- 
serve the presence of extended jointing zones, easily 
detectable in satellite images even up to 200 m res- 
olution. Nevertheless, how well lineaments can be 
detected strongly depends on a number of factors. In 
particular, it depends on the current level of stress 
in the crust. Generally, an enlarged presence of lin- 
eaments indicates that in these regions the crust is 
more permeable, allowing the elevation of fluids and 
gases to the surface. Accumulation of stress deep in 
the crust modifies all afore mentioned processes and 
leads to the variation in the density and orientationn 
of lineaments, previous to a strong earthquake. 
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